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In Abelian theories, one can write for the hard thermal loop equations of motions a local formulation that is 
more economical than the traditional Vlasov formulation and is in explicitly canonical form. I show how this 
formulation can be used for simulating non-equilibrium dynamics in the Abelian Higgs model. 
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1. HARD THERMAL LOOPS 

Understanding the behaviour of quantum fields 
out of equilibrium is important both for cosmol- 
ogy and for heavy ion theory. A particularly in- 
teresting example of an event in the early universe 
in which out-of-equilibrium dynamics played a 
crucial role is the electroweak phase transition, 
since it may explain the observed baryon asym- 
metry in the Universe Although static equi- 
librium properties such as the phase structure of 
the theory have been determined to high accuracy 
with numerical lattice simulations ||^, very little 
is known about the non-equilibrium behaviour. 
The reason is simply that while static correla- 
tors are given by a Euclidean path integral, which 
is well-suited for Monte Carlo simulations, one 
needs to evaluate a Minkowskian path integral to 
obtain any real-time correlator. 

In the case of the electroweak phase transition, 
we can safely assume that the fields are in thermal 
equilibrium before the transition. The occupa- 
tion number of the soft modes (with momentum 
A; -C T) is large, and a classical approximation can 
be used for them [|j . On the other hand, the hard 
modes {k > T) have a small occupation number 
and must be treated as quantum fields. Pertur- 
bation theory works well for the hard modes, and 
we can integrate out them perturbatively, con- 
structing a classical effective theory for the soft 
modes This construction is free of infrared 
problems, since the loop integrals have an infrared 
cutoff, namely the ultraviolet cutoff A of the ef- 
fective theory. This approach has previously been 
used to measure the sphaleron rate of hot SU(2) 
gauge theory [00. It can also be used for simu- 



lating non-equilibrium dynamics of phase transi- 
tions, since the phase of the theory is a property 
of the soft modes only, and the distribution of the 
hard modes does not change in the transition. 

At one-loop level, the resulting classical theory 
is simply the hard thermal loop effective theory. 
In high-tempcrature approximation and neglect- 
ing the IR cutoff in the loop integrals, we end up 
in the expression B 
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(1) 



= m2 + (eV4 + A/3)T2, 



and the integration is taken over the unit sphere 
of velocities v = {l,v), = 1. The time evo- 
lution of the fields is given by the equations of 
motion 
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- 2elm(j)*D''(j),{2) 



(3) 



Because of the derivative in the denominator, the 
gauge field equation of motion is non-local: the 
field interacts with fields on its whole past light 
cone. 

2. LOCAL FORMULATION 

To make numerical simulations practical, one 
needs a local formulation for the theory. For 
non-Abelian theories, two such formulations have 
been proposed: adding a large number of charged 
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point particles Q and introducing an extra field 
W, which satisfies the non-Abelian Vlasov equa- 
tion 1^,0. In practice, both of them lead to a 
S+l-dimensional field theory. However, we will 
use the formulation presented in Ref. [§, where 
it was shown that in the Abelian case one can 
integrate out one of the dimensions, thus obtain- 
ing a 4-|-l-dimensional theory that is completely 
equivalent with the others. It consists of two ex- 
tra fields f{t, X, z) and 9{t, x, z), where z G [0, 1]. 
In the temporal gauge, they satisfy the equations 
of motion 

d^,f{z) = z^v^f + mDZ^^-^VxA, (4) 
d^9{z) = z^V • (V6l - moA^ , (5) 

^TfiuA^ y^^V X / j + J. (6) 

It is now possible to calculate any real-time cor- 
relator at finite temperature. One simply takes 
a large number of initial configurations from the 
thermal ensemble (the explicit form of the Hamil- 
tonian was given in Ref. [^), and evolves each of 
them in time, measuring the correlator of inter- 
est. The average over the initial configurations 
gives the ensemble average of the correlator. In 
principle, that approximates the corresponding 
expectation value in the original quantum theory 
to leading order in gauge coupling constant. In 
practice, one has to worry about some technical 
issues concerning the lattice cutoff that have been 
discussed in Refs. 

3. SIMULATIONS 

To simulate the theory numerically, we dis- 
cretize it by putting it on a lattice [|j . We expand 
/ and B in terms of Legendre polynomials: 

/^"^ = ^ dzZ^^^P2n{z)f{z), 

- f' dzP2n{z)0{z), (7) 
JO 

In addition to equilibrium real-time correlators. 



this formalism can be used to simulate out-of- 
equilibrium phenomena, such as the dynamics of 
the first-order phase transition ||l^ at A <C e^. 
We do this in three steps: 

1. We thermalize the soft modes to T > Tc 
with a Monte Carlo simulation. 

2. We decrease the temperature to slightly be- 
low Tc, and perform a smaller number of 
sweeps and generate the hard modes around 
the soft background. Since we use only lo- 
cal updates, the system does not tunnel to 
the true minimum, i.e., the Higgs phase, but 
remains in the Coulomb phase. 

3. We solve numerically the equations of mo- 
tion using the configuration obtained in 
step 2 as an initial configuration. 

As an example, we show in Fig. |l| how the 
winding number density changes during the phase 
transition. Initially, the system is full of fiuctuat- 
ing strings, but after some time, some bubbles of 
the Higgs phase form and start to expand. The 
bubbles collide and strings are trapped between 
them. This is an example of the Kibble mecha- 
nism [0 . With the technique presented here, we 
can measure the properties of the resulting string 
networks and compare them with estimates based 
on non-gauged theories. 

This approach can also be used to study various 
other properties of the phase transition. One of 
the most interesting quantities is the bubble nu- 
cleation rate. It can be measured in a straightfor- 
ward way by taking a large number of initial con- 
figurations from the supercooled Coulomb phase 
and evolving each of them in time. If we plot the 
number of runs still remaining in the symmetric 
phase after time t, we get an exponentially de- 
creasing curve, whose exponent gives the nucle- 
ation rate. This can be compared with the ana- 
lytical results jl^ to see, how well the analytical 
calculation really works in gauge theories. 
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Figure 1. Time evolution of the winding number 
density during the course of a first-order phase 
transition. 



